Collective modes of a trapped Lieb-Liniger gas: 
a hydrodynamic approach 
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We consider a trapped repulsive one-dimensional (ID) Bose gas at very low temperature. In 
order to study the collective modes of this strongly interacting system, we use a hydrodynamic 
' approach, where the gas is locally described by the Lieb-Liniger model of bosons interacting via a 

repulsive delta potential. Solving the corresponding linearized hydrodynamic equations, we obtain 
the collective modes and concentrate more specifically on the lowest compressional mode. This is 
done by finding models, approaching very closely the exact equation of state of the gas, for which 
the linearized hydrodynamic equations are exactly solvable. Results are in excellent agreement with 
those of the sum rule approach of Menotti and Stringari. 

I. INTRODUCTION 

VO , 

Cold atomic gases of reduced dimensionality are currently being investigated experimentally [1-3]. A typical 
situation is that of a dilute 3D Bose-Einstein condensate (BEC) confined in a very anisotropic trap. If the confinement 
is strong enough in two perpendicular directions, the gas behaves microscopically as a ID Bose gas in the remaining 
perpendicular direction. A very interesting model of such an interacting system is that of Lieb and Liniger [4]: it 
consists of ID bosons interacting via repulsive zero-range interactions of variable strength. In the homogeneous case, 
this model can be solved exactly, whatever the strength of the interactions, by means of Bethe ansatz [4]. 

In the present paper, we propose to study the lowest lying collective modes of the trapped Lieb-Liniger gas by means 
of a hydrodynamic approach [5,6]. The linearized hydrodynamic equations are solved using the recently developed 
1 ■ method of Combescot and Leyronas [5]. The results, already presented in [7], are compared with those of Menotti 
and Stringari [6], which were obtained from a sum rule approach. 
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II. TRAPPED LIEB-LINIGER GAS 



The physical system we consider is a 3D dilute BEC at very low temperature (T ~ 0) confined in an elongated 
cigar-shaped trap. The trap is taken to be harmonic with radial frequency uj±_ = u) x = u> y much greater than the 
axial frequency u) z . At very low temperature, repulsive interactions between atoms are fully characterized by the 3D 
s-wave scattering length a > 0. 

The temperature T and the chemical potential [i are much smaller than the radial quantum of oscillation huj±, so 
that the radial motion is frozen and the gas behaves microscopically as a ID system along the axial direction of the 
trap. Olshanii [8] showed that in that case the corresponding effective ID scattering length a\ is given by a\ = a\ja 
[9] when a± = y^h/muj^ ^> a, which we assume for simplicity. In order to study such a trapped ID Bose gas, we 
assume that it is locally described by the Lieb-Liniger model [4] of hard-core bosons, of mass m, interacting trough a 
repulsive contact potential giS(z), where g\ = 2h 2 /mai [8]. 

In the Lieb-Liniger model, the relevant parameter measuring the strength of the interactions is the dimensionless 
coupling constant 7 = mgi/h 2 ni — 2/n\d\, where n\ is the ID density. It is the only dimensionless intensive 
parameter proportional to the coupling constant g\ [4] . It shows that the gas is weakly interacting 7 <C 1 in the dense 
limit n\a\ 3> 1 and strongly interacting 7 1 in the dilute limit ri\a\ <C 1. 

In the dense limit, the homogeneous gas is a quasi-condensate, i.e. a condensate with a fluctuating phase. Never- 
theless, in a trapped gas, when the temperature is well below the degeneracy temperature [10], the phase is almost 
constant over the whole sample and the gas behaves as a true condensate: it is, for example, well described by a ID 
Gross-Pitaevskii equation [11]. In this limit, we say that the gas is in the "ID mean-field" regime and the chemical 
potential is given by \x = g\ri\ — 2huj±ani (with the convention that fi = when n\ = 0). The low density limit 
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corresponds to the "Tonks-Girardeau" regime, i.e. to a gas of impenetrable bosons. Girardeau [12] showed that in 
this limit there is an exact mapping between the gas of impenetrable bosons and the ID free Fermi gas. In the 
Tonks-Girardeau regime, the chemical potential is given by ji = (hirni) 2 /2m. In the whole range between the ID 
mean- field and the Tonks-Girardeau regime, the equation of state /u(ni) can be obtained from the exact solution of 
the homogeneous Lieb-Liniger model [4] . Mcnotti and Stringari numerically evaluated this equation of state and made 
their result available [6,13]. As we will see, the equation of state is the essential ingredient needed in the hydrodynamic 
approach in order to compute the collective modes. 

In the following, we use reactive hydrodynamics, where dissipation is negligible and thermal effects can be omitted. 
This is valid at low enough temperature. The hydrodynamic equations are: 

^+d z {n lV ) = i) (1) 

where ni(z, t) is the ID density at position z and time t and v(z, t) is the velocity along the z axis. The first equation 
is the continuity equation, the second is the Euler equation at zero temperature. The equation of state for the 
homogeneous gas /x(ni) appears as a needed input in the hydrodynamic equations: it is provided by the Lieb-Liniger 
microscopic theory. 

Here we briefly discuss the validity of the hydrodynamic approach. First, in the ID mean-field regime, a trapped 
gas with a finite number of particles behaves as a true condensate at low enough temperature [10]. It is therefore 
well described by a ID Gross-Pitaevskii equation. As in the 3D case [14], neglecting the quantum pressure (which 
is valid when the number of particle is large), the Gross-Pitaevskii equation can be rewritten in the form of reactive 
hydrodynamics with the equation of state given by the mean-field result = gin\. Second, for the whole range of 

densities (i.e. of coupling constant 7), the low lying excitations of the homogeneous Lieb-Liniger gas are known to be 
phonon-like with a sound velocity depending on the coupling constant [4] . This property allows the description of the 
system in term of an effective Luttinger liquid theory [15]. In the case of the trapped gas, when studying collective 
modes of long wavelength and low energy, the chemical potential is locally given by its value in the homogeneous gas 
at the same density /x[rii(z)]. This allows to extend the hydrodynamic approach initially developed for condensates 
in the mean-field regime [14] to the strongly interacting regime [5,6]. 



III. COLLECTIVE MODES : SOLVABLE MODELS FOR HYDRODYNAMICS 



In this section we use the method recently developed by Combescot and Leyronas [5] in order to solve the linearized 
hydrodynamic equations. The basic idea is to replace the equation of state by a model fJ, mo d(ni) for which the 

linearized hydrodynamic equations can be solved analytically or quasi-analytically. The model has several parameters 
that allows to fit the actual equation of state as closely as possible. Details of the fitting procedure are given in [7]. 
We start by briefly reviewing the method of Combescot and Leyronas in the ID case. 

At equilibrium the particle density n\{z) satisfies: 

H[n\{z)] +muj 2 z z 2 /2 = p, (3) 

where jl is the constant value over the system of the global chemical potential. For small density fluctuations 
5ni(z, t) = ni(z, t) — n\{z) we introduce the departure of the chemical potential from its equilibrium value w(z, t) = 
p(z,t) — Jh = (dfi/dn ( l)Sni(z,t). For a fluctuation occurring at frequency u, the preceding equation gives iujSni = 
d z (n\v) — n\d z v + vd^n®) which, together with Euler equation (2) leads to: 

n\d 2 z w + d z n\d z w + mu 2 (dn° 1 /diJ,)w = (4) 

The equilibrium relation (3) implies (dy,/ dn\)(drii/dz) = —muj 2 z which gives: 

zw" + zL'(z)w' - v 2 L\z)w = (5) 

where v 2 = ui 2 /u) 2 is the squared reduced mode frequency and we defined L(z) = ln[n°(z)] with L'(z) = dL/dz. Let 
us introduce the function v(z) by w(z) — z l v(z), with I — for an even mode with respect to z and I — 1 for an odd 
mode, and restrict to z > 0. With this change of function, equation (5) becomes: 

zv" + [21 + zL'{z)}v' - {v 2 - l)L'(z)v = (6) 
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This equation is invariant under the replacement z — > z/R. We therefore take the cloud radius R as unity in the 
following. The cloud radius is given by ^[n°(0)] = muo 2 R 2 /2. Next we make the change of variable y — z a to obtain: 

d 2 v ^ ^ dL dv v 2 - I dL ^ 
dy 2 dy dy a dy 

where A = 1 + (21 — 1) /a. Finally, we define the reduced density n\(z) = n\(z)/n\(Q) and the normalized local 
chemical potential p,(z) = fj,[m(z)]/ '/j,[ni(0)] — 1 — z 2 . We note that in the mode equation (7), the equation of state 
only enters trough the quantity dL/ dy. 

The simplest model introduced by Combescot and Leyronas is the two parameter a — p model [5] . It consists in 
replacing the quantity dL/dy by the model — p/(l — y) in equation (7), which then reduces to the hypergeometric 
differential equation: 

d 2 v dv v 2 — / 

V{1 V) dy^ + [A - V{P + A)] dy- +P ^T V = ° (8) 

The solution to this equation (with proper boundary conditions) is a polynomial [5]. The corresponding mode 
frequencies are: 

o , oe 21 — 1. 

v 2 =l + -n(n+p + ) (9) 

p a 

where n is a positive integer and I = or 1. The a — p model corresponds explicitely to the equilibrium density 
fii(z) = (1 — z a ) p and to the equation of state ji = 1 — (1 — n^ p ) 2 l a . We can check on equation (9) that the "dipolc 
mode" (n = and I = 1) occurs at frequency lo = u) z in agreement with Kohn's theorem. The particular case a = 2 

corresponds to the power law equation of state p, = n^ v . This is precisely the functional dependence of the equation 
of state in the two limiting regimes that we discussed above: the Tonks-Girardeau regime corresponds to p = 1/2 and 
the ID mean-field regime to p = 1 . In the Tonks-Girardeau case, the mode frequencies are given by u> = ku> Zl where 
k = 2n + 1 is any positive integer (it is convenient to notice that I 2 = I). This result was already obtained in [16,6]. 
In the ID mean field regime, the mode frequencies are given by v 2 = k(l + fc)/2, where k = 2n + I is any positive 
integer, in agreement with [6]. 

Combescot and Leyronas introduced other models which lead to a generalized hypergeometric differential equation 
with quasi-polynomial solutions [5]. Here, we only describe the "3 parameters quasi-polynomial" model. The quantity 
dL/dy is modeled by — (p + piy)/(l — y) depending on the three parameters a, p and p\. It is solution of: 

yi 1 -y)i-^ + (-piy 2 -i A + po)y + A )-r- + (pi» + = ( 10 ) 

dy z dy a 

This differential equation admits solutions that are very rapidly converging series which can safely be truncated above 
some order and that we call quasi-polynomials [5] . The mode frequencies are easily obtained numerically. Details are 
given in [7]. 

From now on, we will concentrate on the lowest compressional mode (n = 1 and 1 = 0). In order to obtain the 
mode frequency, we will use the following procedure: (i) fit the actual equation of state by a model (either a — p or 
3 parameters quasi-polynomial) to obtain the best set of parameters; (ii) calculate the mode frequency using this set 
of parameters (with equation (9) in the case of the a — p model, for example). In the next section, we show how to 
improve on this zero-order result for the mode frequency by using perturbation theory. 

IV. PERTURBATION THEORY 

In order to account for the difference between the actual equation of state and the model equation of state we can 
rewrite the mode equation (7) in the form of a Schrodingcr-like equation and use standard perturbation theory. This 
allows to correct the results obtained for the mode frequency calculated using the model equation of state and to 
come closer to the mode frequency corresponding to the actual equation of state. 

The following change of function ij)(y) = v(y) [y A ni(y)] is made in equation (7) leading to: 



dy 2 + 



2 4 K a 2 ' y 4y 2 
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which is a Schrddinger equation with an effective potential (term between brackets [ ]) and corresponding to zero 
energy (and h 2 /2m = 1). The effective potential can be decomposed into a zero-order term corresponding to the 
model equation of state and a perturbation term corresponding to the difference of effective potential between the 
actual and the model equation of state. The perturbation is treated using standard first order perturbation theory. 
This leads to a first order change of the energy in the Schrddinger equation. In order to keep the energy equal to 
zero, we have to give a compensating variation 8v 2 for the frequency: 

Su 2 = a Jo y . — — !-S 12 

/„%(-<) i/*- 1 * 2 

where SL' = L' — L' mod . This result is proven in [7]. 



V. LOWEST COMPRESSIONAL MODE 



We apply now the previously developed method to study the lowest compressional mode (n = 1 and 1 = 0), which 
physically corresponds to a breathing motion of the gas along the axis. The needed input is the equation of state 
n(ni) of the Lieb-Liniger model. Using the numerical evaluation of this equation of state by Menotti and Stringari, we 
compute the mode frequency as follows: for each value of the dimensionless density at the center of the trap m(0)ai, 
(i) we fit the equation of state (for ni(z) varying between and ni(0)) with either the a — p or the 3 parameters 
quasi-polynomial models; (ii) the zero order mode frequency is then obtained by inserting the value of the best set of 
parameters in the formula giving the mode frequency for the model; (iii) we then compute the first order perturbation 
correction to the mode frequency with equation (12). 

Figure 1 shows the results of four different calculations of the lowest compressional mode using respectively the 
a — p model and the 3 parameters quasi-polynomial model without perturbative correction, and the same models 
corrected by first order perturbation theory. Actually, at this scale, it is hard to see any difference between these 
calculations, except for the a — p model which gives a slightly lower mode frequency. 
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FIG. 1. Lowest compressional mode of the trapped Lieb-Liniger gas. 



The squared mode frequency v 2 = u> 2 /lj 2 is plotted as 
a function of log 10 [ni(0)ai]. The crosses (+) correspond to the a — p model; the stars (*) to the corrected a — p model; the 
circles (o) to the 3 parameters quasi-polynomial model; and the diamonds (o) to the corrected 3 parameters quasi-polynomial 
model. 



In order to compare the various calculations, we plot the same data in a different manner, see figure 2. As discussed 
in [7] we expect the corrected 3 parameters model to be the most precise of our results. We therefore take this 
corrected model as a reference and plot the differences between the mode frequencies calculated within the three 
other approaches (a — p , corrected a — p and 3 parameters quasi-polynomial models) and this reference, as shown in 
figure 2. We can see on this figure that, once corrected, the results of the a — p and the 3 parameters quasi-polynomial 
models agree remarkably well, at the absolute 10~ 3 level for v 2 . As already mentioned, the a. — p model is exact in 
the ID mean-field and in the Tonks- Girardeau limits where it predicts v 2 = 3 and v 2 = 4, respectively. In between 
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these limits, the uncorrected a — p model has some difficulties in predicting precisely the correct value of the mode 
frequency. For a discussion of this issue see [7] . 
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FIG. 2. Lowest compressional mode of the trapped Lieb-Liniger gas. The squared mode frequency v = u> /lu z calculated 
using the corrected 3 parameters quasi-polynomial model is taken as a reference. The difference between this reference and 
three different calculations (a — p, corrected a — p and 3 parameters quasi-polynomial models) are plotted as a function of 
log 10 [ni(0)oi]. The crosses (+) correspond to the a —p model; the stars (*) to the corrected a — p model; and the circles (o) 
to the 3 parameters quasi-polynomial model. 

Next we want to compare our most precise result (the corrected 3 parameters quasi-polynomial model) with the 
results of Menotti and Stringari [6], who calculated the lowest compressional mode frequency using a sum rule 
technique [14], which is basically a variational approach. This method is known to provide an exact upper bound to 
the lowest mode frequency. Both results are shown in figure 3. It is seen that they agree very well over the whole 
range of m(0)ai. 
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FIG. 3. Lowest compressional mode of the trapped Lieb-Liniger gas. The squared mode frequency v 2 = lo 2 /to 2 is plotted as 
a function of log 10 [ni(0)ai]. The diamonds (o) correspond to the corrected 3 parameters quasi-polynomial model, and the full 
line is the result of Menotti and Stringari. 
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VI. DISCUSSION 



In this paper, we used a recently developed method to solve the linearized hydrodynamic equations of the trapped 
Lieb-Liniger gas at T ~ 0. We obtained the collective modes with high precision and especially studied the lowest 
compressional mode which is relevant to current experiments. Moritz et al. [3] recently measured this mode in a 
trapped ID Bose gas in the mean-field regime. They found v 2 ~ 3 when the gas is degenerate, in agreement with our 
calculations. As first emphasized in [6], measuring the lowest compressional mode should allow to reveal the transition 
from the ID mean-field to the Tonks- Girardeau regime, where fermionization is expected. 

More generally, we showed, in the particular case of the trapped Lieb-Liniger gas, how the hydrodynamic approach 
allows to obtain collective modes even in a strongly interacting system. We also checked that the sum rule approach 
used by Menotti and Stringari [6] indeed gives very good results for the trapped ID Bose gas. For a discussion of a 
system where this is not necessarily true see [7]. 



ACKNOWLEDGMENTS 



We are grateful to Sandro Stringari, Chiara Menotti and Dima Gangardt for many stimulating discussions. Labora- 
toire Kastlcr Brosscl and Laboratoire de Physique Statistique are Unites de Recherche dc l'Ecole Normalc Superieure 
(ENS) et des Universites Paris 6 et Paris 7, associees au Centre National de la Recherche Scientifique (CNRS). 



[1] Schreck, F., et al, 2001, Phys. Rev. Lett., 87, 080403. 
[2] Gorlitz, A., et al, 2001, Phys. Rev. Lett., 87 , 130402. 

[3] Moritz, H., Stoferle, T., Kohl, M., and Esslinger, T., 2003, arXiv:cond-mat/0307607. 

[4] Lieb, E.H., and Liniger, W., 1963, Phys. Rev., 130, 1605. 

[5] Combcscot, R., and Leyronas, X., 2002, Phys. Rev. Lett., 89, 190405. 

[6] Menotti, C, and Stringari, S., 2002, Phys. Rev. A, 66, 043610. 

[7] Fuchs, J.N., Leyronas, X., and Combescot, R., 2003, arXiv:cond-mat/0305647, to appear in Phys. Rev. A. 
[8] Olshanii, M., 1998, Phys. Rev. Lett., 81, 938. 

[9] We use a definition of the effective ID scattering length similar to that of Menotti and Stringari [6] . There is a difference 
in sign with the definition of Olshanii [8]. 
[10] Petrov, D.S., Shlyapnikov, G.V., Walraven, J.T.M., 2000, Phys. Rev. Lett., 85, 3745. 
[11] Gross, E.P, 1961, Nuovo Cimento, 20, 454; Pitaevskii, LP., 1961, J. Exp. Theor. Phys., 40, 646. 
[12] Girardeau, M., 1960, J. Math. Phys., 1, 516. 
[13] We thank Chiara Menotti for providing us with her data. 

[14] Dalfovo, F., Giorgini, S., Pitaevskii, L.P., and Stringari, S., 1999, Rev. Mod. Phys., 71, 463. 
[15] Haldane, F.D.M., 1981, Phys. Rev. Lett., 47, 1840. 

[16] Minguzzi, A., Vignolo, P., Chiofalo, M.L., Tosi, M.P, 2001, Phys. Rev. A, 64, 033605. 



G 



